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ABSTRACT 



A method for nonlinear system design, based on the generalized 
Popov stability criterion, is developed. It is shown that, although the 
design method is applicable for all values of q in the Popov theory, it is 
particularly useful when q is non-zero, since it greatly simplifies the 
design effort for this case. All designs are accomplished in conjunction 
with the Nyquist and modified Nyquist loci. Basic to the design proce- 
dure is the development and utilization of the modified frequency-response 
polar loci for lag and lead compensation networks. Three examples, one 
with digital simulation, are included to illustrate the procedure described. 
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I . INTRODUCTION 



The now well-known Popov criterion has provided a new approach 
to the idea of absolute stability in the dynamics of automatic control 
systems. As a result, the large amount of subsequent research has 
contributed much to the development of modern stability theory. 

The importance of stability theory is quite evident when one realizes 
that all the current conventional design techniques for linear, time-invar- 
iant systems are directly or indirectly derived from stability theory. This 
is reflected in the fact that the basic engineering design approach is to 
determine how an unstable system can be altered to make it stable. 

When the system is nonlinear or time-varying, the design problem 
becomes much more complex. However, considerable progress has been 
made in this area in recent years . 

This paper considers systems of the standard form of Figure 1. The 
linear element may contain time delays or distributed parameters. The 
nonlinear element may be time-varying and may contain hysteresis . The 
only restriction on the nonlinear element is that its input-output character- 
istic must lie within a sector bounded by two straight lines passing 
through the origin. There is no loss of generality if one line corresponds 
to the x-axis since this can be accomplished through a simple transfor- 
mation. Hence: 

0 ^ £ k for all t. ( 1 ) 

e (t) 
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Fig. 1. The basic system considered in this paper. 



This class of systems may be treated, in some cases, by phase- 



plane methods or by approximate methods such as the single or dual- 
input describing function. Some exact methods are also historically 
available . 

The solution of the absolute stability problem was first formulated 
in 1944 by Lur'e and Postnikov [1] for the case when u = f(e) is located 
completely in the first and third quadrants (k = °°). Also, in their paper, 
the utilization of a Lyapunov function of the type "a quadratic form plus 
an integral involving the nonlinearity" was first proposed. In further 
works, Lur'e developed a method which permits one to write down a 
finite system of quadratic equations directly from the given differential 
systems. He showed that if this system of equations, which he called 
resolving equations, has real roots, then this is a sufficient condition 
for the absolute stability of the system of Figure 1, 

Yakubovitch [1] further developed this approach by showing that the 
system of resolving equations can be reduced to a single algebraic 
equation. He then proved that the existence of two real roots for this 
equation is a sufficient condition for absolute stability. 

Most of the subsequent research was centered around the application 
of Lyapunov's method for determining system stability. However, up 
until this time, all the exact methods suffered from the disadvantage 
that even the solution of the simple problems became so formidable 
computationally that the methods were generally unusable as engineering 
tools . 
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Then, in 1959, Popov proposed an entirely different approach to the 



problem of absolute stability. He expressed his sufficient conditions 
in terms of the frequency response of the linear portion of the system, 
thus giving his criterion a very simple form that is convenient for graph- 
ical application. In fact, the Popov criterion is so powerful that all the 
results described above and connected with the Lyapunov function, con- 
sisting of a quadratic form plus an integral of the nonlinearity, are 
included in the criterion. 

In the Popov method only the Nyquist diagram of the linear element 
and the inequality constraints of the nonlinear element are needed. 

Hence the criterion is readily applied to higher-order systems and those 
containing time delays or distributed parameters . Since it is basic to 
the design method developed in this paper a brief development of the 
Popov theorem follows. For a proof of the theorem, see [l] . 

In the most simple case, the problem of absolute stability of dynamic 
processes consists of finding when the equilibrium state of the general 
system 

x = Ax + bu 

u = f (e) (2) 

e = c'x 

where x is a vector, A is a square matrix, b is a column vector, c' is a 
row vector, and u is a scalar, is on the whole asymptotically stable in 
the Lyapunov sense. The formulation (2) is a system of differential 
equations differing from a linear homogeneous set with constant 



10 



coefficients only by a nonlinear function whose argument is any one or a 



linear combination of the state variables. 

A, b, and c are real, time -invariant matrices in which some of the 
components may be zero. The nonlinear characteristic, f(e) , is an 
arbitrary, single-valued, piecewise continuous real function, defined 
for all real values of e and satisfying the conditions 

f(0) = 0 (3) 

0 £ ^ k 

e 

which are interpreted geometrically in Figure 2. 

The problem of establishing the absolute stability of the system (2) 
requires the following definitions. A system (2) such that all the roots 
of its linear portion characteristic equation lie in the left half-plane will 
be the so-called principal case of (2). However, a case when some of 
the roots of the characteristic equation are on the imaginary axis and the 
rest are in the left half-plane will be called a particular case of (2). In 
(2), assume that f(e) = he with 0 £ h =£ k , for the principal case and 
with 0 < h ^ k for the particular cases . Then (2) becomes the linear 
system 

k = Ax + hbe 

(4) 

c = c's . 

Now, Popov expressed his results in the frequency domain. In addition 
to the usual open-loop frequency response, G(joj) (see Figure 1), he used 
a modified frequency response, G*(ju>), which is defined as follows 
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Fig . 2 . The input-output characteristic of the nonlinear 
element. 
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Re G*(jw) = Re G(jco) = X 

( 5 ) 

Im G*0'co) = 03 Im GGoj) = Y . 

From linear theory, for the system (4), obtained from the principal 
case of (2), to be stable for any h in the interval 0 ^ h ^ k, it is 
necessary and sufficient that the locus of G(joi) , and hence also the 
locus of G* G 03) , never intersect the real axis from negative infinity up 
to and including the point (-1/k) . Also for the system (4), corresponding 
to particular cases of (2), to be stable with 0 < h ^ k , it is necessary 
and sufficient that the locus of GG 03 ) (or of G*Gw) ) never intersect the 
forbidden zone (the interval of the negative real axis from negative 
infinity to - 1/k) and that the system (4) be stable-in-the-limit , that is, 
that the system (4) be stable for sufficiently small h > 0 (that is, 
h = € > 0 where € is arbitrarily small) , 

Consider now the original problem of establishing the absolute 
stability of the system (2), that is, when f(e) is an arbitrary (not 
necessarily linear) characteristic contained in the sector (0,k) for the 
principal case, or in the sector (e , k) for the particular cases. 
Additional conditions must be imposed on the frequency response in 
order to guarantee stability for the system (2), not only for linear charac- 
teristics f(e) = he , but for any nonlinear characteristics f(e) contained 
in the given sector. Popov stated these additional requirements in the 
following famous theorem [l ] . 

Theorem . For the system (2) to be absolutely stable in the sector 
(0,k) for the principal case, and in the sector (e , k) for the particular 
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cases (where €> 0 is an arbitrarily small number), it is sufficient that 
there exist a finite real number q such that for all co s 0 the following 
inequality is satisfied 

Re(l + jqu>)G(jaO + l/k>0 (6) 

and, additionally for the particular cases, that the conditions for 
stability-in-the-limit be satisfied. 

In order to give the criterion a geometric interpretation let 

G*(jto)=X + jY (7) 

Then 

Re(l + jqo))G(ju>) = ReG(ju>) - qwImG(jw) = X-qY. 

Hence (6) can be written as 

X-qY + 1/k > 0 for all « * 0 . (9) 

But the equation 

X-qY+l/k=0 (10) 

is the equation of a straight line (the Popov line) with slope 1/q which 
passes through the point -1/k on the real axis. Hence it is possible to 
give a geometric interpretation to the Popov theorem. 

For the system (2) to be absolutely stable in the sector (0,k) for the 
principal case, or in the sector (e , k) for the particular cases, it is 
sufficient that there exists in the G*0 cp) plane a straight line, passing 
through the point -1/k on the real axis, such that the modified frequency 
response G*(jo>) lies strictly to the right of it. Also, for the particular 
cases, the conditions for stability-in-the-limit must be satisfied. Some 
examples are shown in Figure (3). 
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Fig. 3. Examples of the application of the Popov criterion. 



Attempts to extend the Popov theorem have met with considerable 
success. Hsu and Meyer [2] have presented a generalized Popov theorem 
which includes the recent extensions . 

Theorem . Consider the basic system of Figure 1 with the linear 
element being output stable. For the system to be both absolutely 
control and output asymptotic for (u/e) € (0/k) , it is sufficient that there 
exists a real number q such that for all real a> ^ 0 and an arbitrarily 
small 6 > 0, the following inequality is satisfied: 

Re ( (1 + jooq)G(jw) ) + 1/k s 6 > 0 (11) 

The restrictions on qand k depending on the nature of the nonlinear 
element, are: 

1. for a single-valued, time -invariant element: 

if 0 < k < 00 , then - »< q < “ 

if k = <® , then 0 ^ q < <® 

2. for a nonlinearity with passive hysteresis: 

0 < k < 00 and - ® < q ^ o 

3. for a nonlinearity with active hysteresis: 

0 < k ^ <=° and 0 ^ q < » 

4. for a general nonlinearity (time -varying , and possibly with 
hysteresis) : 

0 < k ^ 00 and q = 0 

If k = “ in cases (1), (2), or (4) above, then an additional requirement 
must be met. For all t and for every value e B , there must exist a finite 
value u B such that 

I u(t) | ^ u, < » if | e(t) | ^ e„ (12) 

It can be noted immediately that the term "absolutely stable" has 
been replaced by "absolutely control and output asymptotic". This is 



16 



mathematically the least restrictive practical stability definition. It 
simply means that the control signal or the output signal eventually go 
to zero. 

A valuable engineering tool provided by this theorem is the trade-off 
available between the requirements on the linear element and those on 
the nonlinear element. For example, if the nonlinearity is single-valued 
and time invariant, a non-zero q may be chosen, thus providing a greater 
range for G(joo) . 

The Popov criterion provides another valuable asset to the engineer. 
A measure of the rate of system damping is easily incorporated in the 
basic theorem. If the linear element of the system of Figure 1 is output 
stable of degree a , this simply means that, for a > 0 , the system 
output, in response to either an initial condition or an impulse, tends to 

— Q/ t 

zero faster than the function e . For a<0, it means that the output 
OL t 

multiplied by e ultimately goes to zero. 

Now it can be said that if the system of Figure 1 is control asympto- 
tic of degree a. , there exists a real number a such that for every set of 
initial conditions 

CO 

' (e at u(t) ) 3 dt < « • (13) 

0 

Also, it is output asymptotic of degree a if there exists a real number a 

such that for every set of initial conditions 
00 

J (e at e(t) ) s dt < « . (14) 
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In order for the Popov inequality to reflect a measure of damping, the 
linear element must be output stable of degree a , and G(jw-a) must be 
substituted for G(jcc). Hence, for the system of Figure 1 to be control 
and output asymptotic of degree a . , the following inequality must be 
satisfied: 

Re ( (l+jqw)G(jw-a) ) + 1/k s 6 > 0 (15) 

Methods of system design based on the Popov criterion have been 
developed by several authors. See, for example, [2] , [3] , and [4] . 
These methods are all relatively straightforward and easily applied for 
the case q = 0, beacuse of the invariance of the Popov line with frequency 
on the Nyquist plot. However, for the case q / 0, the design of lag-lead 
compensating networks involves considerable design effort, much of it 
devoted to trial -and -error computation. 

A great deal of this labor is eliminated in the new design method 
proposed here for the large group of systems in which design requirements 
dictate using the case q / 0. In order to retain the unparalleled utility 
of the polar-plot representation, the lag-lead compensating networks 
are designed on the G*(jo.') plane. 
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II. DESIGN METHOD 



A sample design objective is first illustrated in terms of a particular 
system. Then the general design method is developed. 

Assume that the linear element of Figure 1 is output stable and has 
the Nyquist plot shown in Figure 4. Suppose also that the system design 
requires the nonlinear element to be a saturating amplifier with the input- 
output characteristic of Figure 5. 

The generalized Popov criterion for a general nonlinearity (q=0) gives 
a maximum value of k = 1.67 for the upper bound on the Popov sector to 
guarantee control and output asymptoticity of the system. This is shown 
in Figure 4. However, since the particular nonlinearity of this system 
falls under Case 1, a non-zero q may be chosen to obtain a less con- 
servative result. In order to do this, the simple transformation is made 
from the G(jU') plane to the G * (j 4*-' ) plane by multiplying the ordinates of 
the G(j£*') curve by their respective frequency values. The result is the 
modified frequency response, G*(jaO. Assume that G*(jU.’) has the form 
shown in Figure 6. Drawing a tangent line to the G*(j&>) curve (see 
Figure 6) at the intersection of the G*(jk') curve and the negative real 
axis gives a value for q of 0.1 and a value for k of 2.0. 

The system nonlinearity (see Figure 5), however, requires a Popov 
sector with an upper bound of at least 9.0 in order to guarantee a stable 
system with the chosen amplifier. Hence the system must be compen- 
sated. Excessive compensation and a too conservative solution can be 
avoided by choosing a non-zero q. The design of lag or lead cascade 
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Fig. 4. Linear system polar plot. Fig. 5. Saturating amplifier input-output 

characteristic. 




Fig . 6 . Polar plot of the linear element modified frequency 
response . 
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compensation of the G* plane for systems of this type can be accomplished 



by the method described below with an appreciable decrease in labor over 
previous methods. The compensator is to be inserted in the output of the 
nonlinear element in the system of Figure 1. Before proceeding with the 
compensator design, several preliminary relationships which are instru- 
mental to the design method are given. 

One of the basic difficulties of designing in the G*(jO>) domain is the 
fact that 

[G(ju>)G c (ju>)] * ^G*(jto)G c *0w) (16) 

where G 0 (ja>) is the transfer function of the compensator. In order to 
show later how this difficulty may be minimized, the exact relationship 
is derived. 

G(ju.') = ReG(jw) + jlmGCiu?) 

G c (ju>) = ReGcOw) + jImG c (jO)) 

(17) 

G*(joO = ReG(jco) + ju>ImG(ju') v 

G 0 *(ja') = ReGcOto) + jOJlmG c (jw) 

Multiplying G(ja))G c (j'a>) , 

G(ja')G c (ju>) = ReG(ju>)ReG c (j6e) - ImG(ja.')ImG c (ju.') 

(18) 

+ j [ReG(jaj)lmG c (ja>) + ReG 0 (j w)ImG (j to) ] 

[G(ja')G c (ja))3 * = ReG(ja))ReG c (ja') - ImG(j to)ImG 0 Ga)) 

(19) 

+ ja'[ReG(jto)ImG c (ja.') + ReG c (jco)ImG(ja})] . 

Multiplying G* (j U>)G 0 *(j a.') , 

G*(jw)G c *0a') = ReG(jU')ReG c (jw) - cc 2 ImG(ja3)ImG c (ja)) 

( 20 ) 

+ j w[ReG(ju.')ImG c (ja)) + ReG c (j 6o)lmG(j 0 -') ] . 
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Comparing (19) and (20), 



[G(jto)G c (ju>)3 * = G*(Jw)G c *Gw) + (o3 2 -l)ImG(j6L0lmG o 0'^) . (21) 

The frequency-response plots of the compensation networks are 
essential to the design procedure. The transfer function of a lag compen- 
sation network is given by 



G c (s) = \ ~~ 2 ~ where t 1 >t 2 (22) 

1 i S T i 

and that of a lead network is given by 

G c (s) = * + STl where n > T 3 . (23) 

1 + st 2 

The Nyquist plots of the two networks are shown in Figure 7. 

The geometric configuration of the modified frequency responses, 
G c *(ju>)/ is now derived for the lag filter only. The lead filter gives the 
same result only with t 1 and t 2 reversed. For the lag filter: 



G c *(jw) 



1 + a.’ 3 % t p 

1 + u> 2 Ti 2 



u> 2 (t 2 - tQ 

1 + 0) 2 Ti 2 



• x = 1 + h) 8 TiT 8 . = <£ (T a “ T X ) 

1 + a« 2 r, 2 ' 1 + u> 2 Ti 2 

Solving for co 2 in the X-expression: 

(1 + W 2 Ti 2 )X = 1 + t*> 2 Tj. t 2 

X + w 2 t-l 2 X = 1 + a> 2 T a T 2 

aj 8 r x a X - a' 2 Ti t 2 = 1 - X 



0! 



2 _ 



1 ~ X 
T: S X - t x r 2 



(24) 



(25) 



(26) 
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Fig. 7. Lag and lead compensation network polar plots. 



Substituting into the Y-expression: 



1 - X 



Y = 



1"l X - T X T; 



(t S - T 1 ) 



1 + 



1 - X 



Ti X - T i T 2 



(1-X) (T a - T,) 

Ti 8 X - Ti T 2 

Ti a X - T x Ta + (l-Xjr^ 

Ti S X - Tj. T 2 

(1 -X) (Ta - n ) 

Ti 2 X - t x t 2 - r x 2 X + r x 2 

(1 “ X) (T 2 — T x ) 

Ti (T X - T 2 ) 

(X-l) (t 2 - Ti ) 

Tl (Ts-Tj 

Therefore , 

Y=l/r 1 (X-l). (2 7) 

Similarly for the lead filter, 

Y = l/r 2 (X-l) , (28) 

which is, of course, the equation of a straight line. This very conven- 
ient result is instrumental to the design method. In order to compute the 
end points of the line (24) is simply evaluated at u.'= 0 and oo= 00 . The 
results for both filters are shown in Figure 8. The frequencies u> = 1/t x 
and co = 1 /t 2 occur at the geometric means of the respective polar plots. 

The ideas developed so far can now be applied to the system with 
the modified frequency response shown in Figure 6. Suppose we 
decide to compensate this system with a lag compensator to force the 
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compensation networks. 



compensated locus, [G(ju>)G 0 (j a?) ] *, to intersect the negative real 
axis at a point £ - 0.1 and at a certain crossover frequency, w 0 . 

This would give a Popov sector of at least (0, 10) which would guar- 
antee a stable system with the required amplifier, provided that a 
Popov line of non-zero q can be drawn through the point (-0. 1,0). 

The procedure to accomplish this alteration of the G*(jw) locus 
is first outlined briefly, then in more detail. The desired negative 
real axis crossover point for the compensated locus is designated as 
P 0 . Then the l/G*Goc) and the G 0 *(jw)/P 0 loci are constructed so 
as to intersect at the point Pi where the value of the 1/G*G w) 
locus corresponds to the chosen crossover frequency, w c , and the 
parameters, r x and r 2 , are then computed so that the frequencies of 
the two curves are equal at the intersection point. This construction 
designs the compensator since it guarantees that 

l/G*(jto) = G c *0w)/P o at ic= oJc 
or 

G*(jU>)G c *Gw) = Pc at = u> 0 

where G*(j oo)G c *(jw) and [G c Gio)GGol') ] * are related by (21). This 
relationship is discussed in detail later. 

To show how easily this is accomplished in practice, the example 
of this section is used. Let P 0 = -0.1 and a' c = 1.0. The first step 
would be to construct the l/G*Gw) curve. However, this is unnecessary 
since only one point of the l/G*Gw) curve is actually used, the one 
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corresponding to oj c = 1.0. The coordinates of this point, P lf are easily 
computed and plotted by inverting the complex number read off the G*(ju.') 
locus for a! 0 = 1.0. 

The next step is to construct the G 0 *(jco)/P 0 Pl°t by simply rotating 
the G c *(jto) plot by 180° and placing the oo= 0 point of the compensator 
at 1/P„. The slope of this straight-line plot is adjusted so that the line 
passes through the point P x . Thus P x corresponds to the intersection of 
the l/G*(co) and the G£(ju!)/P 0 loci, and the frequency at this point 
must be the same, i.e., w= U) c » on both curves. Refer to Figure 9. 

Equation (27) gives the equation of the lag-compensator plot. By 
comparing this with the standard point-slope form of a straight line, it 
can be seen that the slope is 1/t 1 . Hence can be determined by 
simply computing the slope of the line graphically. Then t 2 is 
calculated from either part of (25) using w>= 6l' c . The coordinates, X 
and Y, are simply given by the real and imaginary parts of P x /P t . 

This is a very easy calculation since the value of the squared parameter 
(in this case, r x ) is always substituted numerically, whether the 
filter is lag or lead. Using w= OJ c in the calculation of t 2 guarantees 
that the frequency of the G c *(jw)/P 0 function at P x is the same as that 
of l/G*(jto) at P x . 

To avoid a sign error when computing P x /P 0 , it should be noted that 
X and Y are always positive for a lead compensator and X is positive and 
Y is negative for a lag compensator. It should also be noted that there 
are slope limitations on the straight-line compensator plots. Infinite 
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Fig „ 9. Lag network compensation design. 



slope corresponds to the pole and zero which cancel each other while 
zero slope corresponds to infinite pole-zero separation. 

The frequency response curve of interest in the system design is the 
[G(ja)G c (jto)] * curve. The plot given by the compensation method 
explained above gives the G*(joo)G 0 *(jto) curve. Equation (21) relates 
these two quantities. Note that since U) 0 was chosen to equal 1.0 in the 
above description, the two curves cross the negative real axis at the 
same point. This is not the case in general. 

Assume that the two resulting curves are as shown in Figure 10. 

Here a non-zero q line can be drawn through the point (-0.1,0) and the 
system design is completed. In an actual design problem there may be 
other performance specifications, such as minimum bandwidth, which 
must be met before the design is complete. These additional requirements 
will influence the choice of to 0 and P„ 

In order to generalize the method to any u.' c , (21) is repeated for 
reference . 



[G(jO))G c (jw)] * = G*(jU')G c *(jto) + (U!M)ImG(ju>)ImG c (ju>). (21) 
For simplicity, let (to 3 -l)lmG(jw)ImG 0 (jw) = T. It can be noted 
immediately that 

1 



T = 0 when u = — | and » since ImG 0 (ju>) = 0 (29) 

value for which ImG(jU>) = 0 . 

For the compensated system to be worsened by the effect of T, (i.e, 
the [G(jw)G c (jw)] * curve be more negative than the G*(ju))G c *(jw) curve 
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compensated system. 



31 



at the negative real axis crossover), the term T must be negative for 

W= W 0 . 

For a lag network, since ImG c (jto) equals zero or a negative value for 
all «, this can occur only under the following two conditions: 
if 60 c < 1 and ImG(jU3 c ) < 0 

or if C0 c > 1 and ImG(jco 0 ) >0 . (30) 

For a lead network, since ImG 0 (jU>) equals zero or a positive value for 
all o>, T can be negative only under the following conditions: 
if 0 > c < 1 and ImG(jU! c ) > 0 

orifo> c > 1 and ImG(jC0 c ) <0 . (31) 

Using the methods outlined above, a rapid design procedure can now 
be stated. Since the subject of compensation is treated in most control 
system design textbooks, no attempt will be made in this paper to sum- 
marize the various helpful hints given by other authors for compensation 
design of linear systems using lag and lead networks. Since the design 
method outlined here involves reshaping a linear-system frequency- 
response curve, it is assumed that the designer has a feel for the steady- 
state and transient effects of compensation on linear system performance 
and is cognizant of the advantages and disadvantages of lag and lead 
compensators. See, for example, (5). With this in mind, the design 
steps may be listed as follows: 

1. Plot the polar curves for G(jco) and G*(jcc) as co varies from 0 
to “ . The G*(ju>) plot is easily made by multiplying the 
ordinates of the G(joo) locus by their respective frequencies. 
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2. From the required design specifications and system compo- 
nents choose the type of compensation (lag or lead) and W c . 

3. Use (30) or (31) and the design specifications for the initial 
choice of P„ . 

4. Design the compensator. 

5. Sketch the G*(jw)G 0 *(j4c) locus by evaluating the compensator 
of Step 4 for a minimum number of frequencies. Sketch, if 
necessary, the [G(ja.')G 0 (ju')] * locus using (21). 

6. Evaluate the design in view of the specified performance 
criteria . 

If multiple sections are required the process is simply repeated 
using the [G(ja')G 0 (ja')] * curve in Step 1. Since the design method is 
virtually identical for lead and lag compensators, a detailed discussion 
of lead network design is left to the examples section. It should also be 
noted that the design method introduced above is quite flexible. No 
additional complexity is introduced when a P 0 is desired off the negative 
real axis. 1/P 0 is simply designated at the desired point and the 
G c *(ju.') line is rotated the required number of degrees. 
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Ill, EXAMPLES 



To illustrate how the preceding results might be used, the following 
examples are given. 

Example 1 . It is desired to compensate the system shown in Figure 
11 so that the output signal will be asymptotically stable of degree 0.1. 

To simply guarantee output asymptoticity with the specified type of 
nonlinearity, the minimum Popov sector must be [0,5] . However, for 
the system to be asymptotically stable of degree 0.1, the Popov line 
(for some q) must intersect the negative real axis at a value £-0.2 and 
must not be closer than 0.1 horizontal units to the G*(ju;) curve at the 
closest point . 

The G(jto) and G*(jco) curves are shown in Figures 12 and 13. A lag- 
type compensator and an w c = 1.0 are selected. To avoid the possibility 
of a borderline design, a P 0 of -0.08 is used, as opposed to 0.1. The 
graphical design of the compensator is shown in Figure 14. The 1/P 0 
point is at 1/-0 . 08 = -12 . 5 . To locate point Pj = 1/G*(j 1) , note on 
Figure 13 that G*(jl) = -0.4-j0.2 from which 1/G*(jl) = -2+jl as shown. 
The parameter values, t x and t 3 , are computed using the slope of the 
G c *(jtu')/P 0 plot (see Figure 14) and Equation (25). The resulting compen- 
sator is given by 



Gc(s) 



1 + 1.59s 
1 + 10.5s 



(32) 



Using a minimum number of points the G*(joo)G c *(jo)) locus (see 
Figure 15) is easily sketched. Then, using (21), with (30) as a check, 
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Fig. 12. G(jw) locus for Example 1. 
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Fig. 13. G*(jw) locus for Example 1. 
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Fig. 14. Lag compensation design plot for Example 1 
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the [G(jw)G c (jU)) ] * locus can be quickly determined. In (21), the 
(a> s -l)ImG(j6o)ImG c (j6o) term can be evaluated by simply picking the values 
of ImG(jo>) off the Nyquist plot of G(ja>) and using the expression for Y 
in (25) divided by to for Im G c (jto). 

In this problem, the G*(jto)G c *(jo>) and the [G(jto)G c (jto) ] * loci have 
the same negative real axis crossover at to 0 = 1 and are virtually identi- 
cal at other points. The Popov line is as shown in Figure 15 guaranteeing 
control and output asymptoticity of degree 0.12. 

Example 2 . Consider the nonlinear control system with time delay 
shown in Figure 16. Hsu and Meyer [2] studied this system and obtained 
the following results: 

1. For a class 2 or class 4 nonlinearity (q = 0) , the Popov criter- 
ion is satisfied for 0 <k s 1,85 (see Figure 17). 

2. For a class 1 or class 3 nonlinearity (q ^ 0), the Popov criter- 
ion is satisfied for 0 < k < 1.98 (see Figure 17). 

Now it is desired to compensate this system to permit a Popov sector of 

[0,5] for a class 1 or class 3 nonlinearity. 

The G*(jCo) locus is shown in Figure 17. Lag compensation is chosen 
with an to 0 = 1.8. Equation (30) quickly shows that if a P 0 is chosen 
which gives a minimum Popov sector, the initial design will not produce 
an unstable system but will actually give a greater degree of stability. 
Hence a P 0 of -0.2 is used. Designing the compensator as before (see 
Figure 18) gives the result 



G 0 (s) 



1 + 1.27s 
1 + 4.7s 



(33) " 
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Fig. 17. G*(joo) locus for Example 2. 
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Therefore: 



1. For a class 2 or class 4 nonlinearity (q = 0) , the Popov criter- 
ion is satisfied for 0 <k ^ 4,45 (see Figure 19). This example 
illustrates, as a sidelight, how easily the design method may 
also be applied to the case when q = 0, although the advan- 
tages are not quite as great as when q ^0. 

2. For a class 1 or class 3 nonlinearity (q / 0) , the Popov criter- 
ion is satisfied for 0 < k s 5 (see Figure 19.) 

The permissible Popov sector will actually be slightly greater than above 

as indicated by (30). If exact sector limitations are desired, these are 

quickly obtained by using (21) and (30) to obtain the [GO 6d)G c (jto) ] * 

locus . 

Example 3 . Consider the actuator or indirect control system of 
Figure 20. Hsu and Meyer [2] also studied this sytem using the gener- 
alized Popov theorem. Their intent was to determine the critical feed- 
back gain F = F c such that for feedback gains F ^ F c , the output signal 
y(t) would asymptotically approach a steady-state value dependent on 
the constant input value, r. 

The modified block diagram of this system is shown in Figure 21. 
This zero-input model is possible because the constant input, r, can 
be combined with the constant-initial-condition-response term with no 
loss of generality . 

Using the modified frequency response function, 



G*(jw) = G 2 *(jo>) - jF 



( 34 ) 



where 
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Fig. 20. System for Example 3. 
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For asymptotic stability, the Popov criterion with k = “ then requires 
that 

ReG*(j U') _ qIniG*(j u % ) = ReG 2 *(juO-q(ImG 2 *(jw)-F) > 0. 

G 2 *(ja') is plotted in Figure 22, Since the ImG 2 *(ju.') has a maximum 
value of approximately 0.006, the Popov criterion is satisfied with some 
q >0 if F > 0.006 and the system is then control and output asymptotic. 
Hence e(“) - 0 and 



y(°°) = 



G i ( o ) 



1 



(36) 



F + G x (0) F + 1 

for this system. Therefore the range of F for which the system is 
control and output asymptotic is extremely important since F directly 
affects the steady-state error. 

Digital simulation of the system, with F = 0.006, provided results 
agreeing with the theoretical values. The response of the simulated 
system to an input value of r = i . 0 is shown in Figure 23 . 

Assume now that design requirements specify a steady-state error 
which requires the system of Figure 20 to be asymptotically stable for 
F = 0.004. Hence the system must be compensated so that the maximum 
value of the ImG 2 *(jo.') is less than 0.004. For comparison purposes, 
the unit step response of the uncompensated system with F = 0.004 is 
shown in Figure 24. Lead compensation, u! c ~ 8.8, and P 0 = - .017 are 
chosen. The compensator design is shown in Figure 25 with the result 



G c (s) 



1 + 0 . Is 
1 + 0.079s 



(37) 
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Fig. 23. Response of the system of Figure 20 to a unit step 
input. F = 0.006. 



50 



y(t) 




0.6- 



0.4- 



0 . 2 - 



t(sec) 

1 1 1 1 -—— 1 — 

10 20 30 40 50 

Fig. 24. Response of the system of Figure 20 to a unit step 
input. F = 0 . 004 . 
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Fig. 25. Lead compensation design plot 
for Example 3 . 
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The G 0 *(ja>)G 2 *(jaj) locus is shown in Figure 26. 

It should be noted that, in this case, the T term of (29) is of no 
consequence since it is only real and will have no effect on the maximum 
imaginary value of the compensated curve. Hence 

maxim [G c (j w)G 2 (jt»)) ] * = maxim G c *(jti;)G 2 *(ja>) (38) 

and the design is completed. The response of the compensated system 
with F = 0 . 004 and r = i . 0 is shown in Figure 27 . 
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Fig. 26. G 2 *(ju))G c *(ja;) locus for Example 3. 
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Fig. 27. Response of the compensated system of Example 3 
to a unit step input. F = 0.004. 
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IV. CONCLUSION 



A design method, based on the Popov stability criterion, has been 
developed for the large class of nonlinear systems which can be modeled 
in the configuration of Figure 1. The design procedure uses the modified 
frequency-response polar loci for the linear portion of the system. This 
permits the engineer to utilize the clear physical picture inherent in the 
Nyquist plot. 

As a major step in the development of the design procedure, it was 
shown that the modified frequency-response polar plots for lag and lead 
compensation networks are straight lines. This provided the basis for 
the development of a rapid, straightforward design method which could 
be quite useful in engineering practice. 
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